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ABSTRACT 

Warped accretion discs in Active Galactic Nuclei (AGN) exert a torque on the black 
hole that tends to align the rotation axis with the angular momentum of the outer 
disc. We compute the magnitude of this torque by solving numerically for the steady 
state shape of the warped disc, and verify that the analytic solution of Scheuer & Feiler 
(1996) provides an excellent approximation. We generalise these results for discs with 
strong warps and arbitrary surface density profiles, and calculate the timescale on 
which the black hole becomes aligned with the angular momentum in the outer disc. 
For massive black holes and accretion rates of the order of the Eddington limit the 
alignment timescale is always short ( < 10 6 yr), so that jets accelerated from the inner 
disc region provide a prompt tracer of the angular momentum of gas at large radii in 
the disc. Longer timescales are predicted for low luminosity systems, depending on the 
degree of anisotropy in the disc's hydrodynamic response to shear and warp, and for 
the final decay of modest warps at large radii in the disc that are potentially observable 
via VLBI. We discuss the implications of this for the inferred accretion history of those 
Active Galactic Nuclei whose jet directions appear to be stable over long timescales. 
The large energy deposition rate at modest disc radii during rapid realignment episodes 
should make such objects transiently bright at optical and infrared wavelengths. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - galax- 
ies: jets - galaxies: nuclei - hydrodynamics 



1 INTRODUCTION 

The angular momentum of gas accreting onto massive black 
holes in Active Galactic Nuclei (AGN) is likely to change 
with time, either as a result of mergers tunneling fresh gas 
towards the nucleus (Hernquist & Mihos 1995), or because 
of radiation or disc wind driven warping instabilities in an 
existing accretion flow (Pringle 1996, 1997; Maloney, Begel- 
man & Pringle 1996; Maloney, Begelman & Nowak 1998). 
In either case, the angular momentum of gas at large ra- 
dius in the disc will be misaligned with the rotation axis of 
the black hole, while at small radius the combined action of 
viscosity and differential precession induced by the Lense- 
Thirring effect (Lense & Thirring 1918) leads to alignment 
of the disc and hole angular momenta (Bardeen & Petter- 
son 1975; Kumar & Pringle 1985). Achieving this alignment 
requires that the black hole exert a torque on the disc gas, 
and implies that an equal and opposite torque act on the 
black hole itself. Over time, this causes a change in the spin 
axis of the hole towards alignment with the large angular 
momentum reservoir provided by the disc at large radius. 

Whether the timescale for alignment is long or short 
compared to the lifetime of an AGN has implications for our 



understanding of the accretion history of these objects. If the 
black hole is spinning, then jets, regardless of whether they 
derive power from the aligned inner disc or directly from 
the hole itself (Blandford & Znajek 1977; see also Ghosh & 
Abramowicz 1997; Livio, Ogilvie & Pringle 1999), trace the 
spin of the black hole. If the alignment timescale is long, then 
we expect that the jet direction reflects the initial spin of the 
black hole established during (or prior) to the formation of 
the host galaxy (Rees 1978). Moreover the jet direction is 
expected to be stable over time, irrespective of variations in 
the angular momentum of accreting gas. Conversely, if the 
timescale is short then jets define the angular momentum 
of inflowing material, and the observation of systems where 
they appear to be stable over long time periods of 10 7 - 
10 8 yr (Alexander & Leahy 1987; Liu, Pooley & Riley 1992; 
Scheuer 1995) implies a constancy in the average angular 
momentum of accreting gas. 

The rate of realignment was computed by Scheuer & 
Feiler (1996), who found an approximate analytic solution 
to the equations governing the evolution of a warp. The 
timescale depends on two poorly known viscosities, v\ and 
vz, which correspond to the (R,(f>) and (R,z) components 
of the shear as defined by Papaloizou & Pringle (1983). If 
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V\ ~ V2 (i.e. the timescale for warp evolution is compara- 
ble to that for the evolution of the surface density) then 
the formula of Scheuer & Feiler (1996) is similar to that 
previously assumed by Rees (1978). Subsequently, Natara- 
jan & Pringle (1998) revisited the problem, and pointed out 
that for the parameters of AGN discs the assumption that 
V\ ~ U2 is at variance with the results of detailed analyses of 
the hydrodynamics of viscous discs, which suggest instead 
that V2 S> v\. Recalculating the alignment timescale with 
this modification they found that it was short — much less 
than any reasonable estimate of AGN lifetime — and used 
this to argue that the spin of black holes in AGN, and the 
direction of jets produced in the inner regions of the accre- 
tion flow, rapidly adjust to trace the angular momentum of 
gas at large radius. 

The analysis in Natarajan & Pringle (1998) employed 
the same torque formula as Scheuer & Feiler (1996), which 
was derived under the restrictive conditions of small ampli- 
tude warps and a disc viscosity that is constant with radius 
(or, equivalently, that in a steady-state the disc surface den- 
sity is independent of radius) . Since AGN may well harbour 
strongly warped discs, and almost certainly have more com- 
plex surface density profiles, our goal in this paper is to gen- 
eralise their results. Section 2 below sets out the equations 
for the evolution of a warped disc, which we solve numeri- 
cally in Section 3 for comparison with the approximate an- 
alytic solution. The resulting torque is calculated in Section 
4, and the alignment timescale derived for a simple AGN 
disc model in Section 5. Section 6 discusses the implications 
of our results. 



2 WARPED DISC EQUATIONS 
2.1 Governing equations 

The evolution of a thin, warped accretion disc in a Keple- 
rian potential was studied by Papaloizou & Pringle (1983). 
For a disc with a surface density profile E(_R, t), with angular 
momentum at radius R parallel to the unit vector I, the evo- 
lution can be expressed most compactly in terms of an equa- 
tion for the angular momentum density L = (GMR) 1 ' 2 SZ. 
Adopting the simplest description (Pringle 1992; see also 
Papaloizou & Pringle 1983), the time evolution is given by, 



to be general power laws with radius, 
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Here v\ and V2 are viscosities acting on the (R, cf>) and (R, z) 
components of the shear, as defined by Papaloizou & Pringle 
(1983). Neither is equal to the shear viscosity 'u' that enters 
into the Navier-Stokes equations, and the general relation 
between these three quantities is extremely complex (Ogilvie 
1999). For the parameters of AGN discs, the simplest as- 
sumptions suggest that V2 S> v\ (Natarajan & Pringle 1998), 
though for most of this paper we will allow both viscosities 
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For thin, vertically averaged accretion disc models, the shear 
viscosity v is a power-law in radius in a given opacity regime 
(e.g. Frank, King & Raine 1992). Later, the index (3 will 
be chosen to correspond to the scaling expected at radii 
in the disc where the bulk of the realignment torque acts. 
However, for a general warped disc, the radial run of both 
v\ and V2 will also depend on the shape of the warp (Ogilvie 
1999). This effect is not taken into account in this paper. 
Our assumption that v\ and V2 are simple power-laws in 
radius, independent of warp shape, is therefore strictly valid 
in the limit of small warps, and is an approximation for 
larger amplitude warps. 

The range of behaviour allowed by equation (|l|) is 
considerably more diverse than simple diffusion of L. 
Nonetheless, we will make frequent use of the characteristic 
timescales for diffusion of surface density and warp, which 
we define conventionally as, 
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The final term on the right hand side of equation ([l]) repre- 
sents the effect of Lense-Thirring precession on the disc. For 
a black hole with mass M, we take uj p = u P (0, 0, 1), where 



o p = 2ac ^ — —\ 



(4) 



where a is the dimensionless angular momentum of the black 
hole. 

Several assumptions and simplifications are necessary 
to derive this equation, and these are set out in detail in 
Papaloizou & Pringle (1983) (see also Pringle 1999; Ogilvie 
1999). Particularly important is the boundary between the 
diffusive behaviour described by equation (jjj) and wave-like 
evolution, which occurs when the usual Shakura & Sunyaev 
(1973) viscosity parameter a ~ H/R, where H is the disc 
scale height (Papaloizou & Lin 1995) . The observational con- 
straints on a in AGN discs are weak (Siemiginowska & Cz- 
erny (1989) and Siemiginowska, Czerny & Kostyunin (1996) 
suggest a = 1CP 1 — 1CP 2 ), but in the ionized parts of the 
disc it is plausible to assume that a ~ 0.1, as in better stud- 
ied disc systems (eg. Cannizzo 1993). In any case, theoreti- 
cal estimates for H/R at the relevant radii (Collin-Souffrin 
& Dumont 1990) suggest that in AGN a > H/R, so that 
we are safely in the diffusive regime. In this respect AGN 
discs are very different from protostellar discs, for example, 
where wave-like behaviour is likely to be important (e.g. Lar- 
wood et al. 1996). We also note that although Papaloizou 
& Pringle (1983) considered only small amplitude warps de- 
scribable using linear perturbation theory, Ogilvie (1999) 
has shown that strongly warped discs can also be described 
using a similar approach. 

2.2 The Scheuer & Feiler analytic solution 

Scheuer & Feiler (1996) obtained an approximate steady 
state analytic solution to equation (Q) for the case where 
(3 = 0, in the limit where the angle of misalignment of the 
disc and the hole spin was small. In this case, the surface 
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density obeys the usual relation for a planar disc with ac- 
cretion rate M and zero torque inner boundary condition 
imposed at Ri n , 




The shape of the disc is given by / = (l x ,l y ,lz), where 
l x = Ke~ cos <f> 

l y = Ke~^ sin <j> (6) 

and (f> — 2(tjp/^2-R) 1//2 - The constant K depends on the disc 
inclination i at large radius, at large radius K = sini. 



2.3 Numerical solutions 

The above analytic solution applies only for j3 = 0, and as 
derived assumes that the disc warp is of small amplitude. 
Relaxing these assumptions requires a steady state numeri- 
cal solution of equation (|l]). 

For /3 < 1/2, steady state solutions can be calculated 
efficiently by solving directly the ordinary differential equa- 
tions (ODEs) obtained by setting dT,/dt = in equation (|l|). 
We give a full description of our methods in the Appendix, 
but briefly we have found it simplest to use equations ex- 
pressed in terms of the surface density E(-R) and unit tilt 
vector I — L/|L|, which can readily be solved iteratively us- 
ing a finite-difference technique described by Pereyra (1979). 
We use the Numerical Algorithms Group's implementation, 
and impose boundary conditions of E(iZi„) = 0, E(i? out ) = 
Bout, t(R = Rin) = (0,0,1), l(R = Rout) = L,t. These 
boundary conditions assume that the inner disc is aligned 
with the spin axis of the hole, we therefore have to check 
post facto that this is indeed a valid assumption for each 
model calculated. 

For j3 > 1/2, we have been unable to obtain a con- 
verged numerical solution to the steady state ODEs using 
these methods. In this regime, we instead evolve the time 
dependent equation (Jj]) from an arbitrary initial condition 
(of a flat, uniformly tilted disc) until a steady state is ob- 
tained, using the numerical method described by Pringle 
(1992). This is straightforward but computationally expen- 
sive, since many viscous times at the outer edge of the disc 
are required to obtain a steady state solution. Defining the 
viscous time of the disc by t Vl = Rout/ u i, we nn d that more 
than 10 viscous times of evolution are necessary. Moreover, 
for either method capturing the torque accurately requires 
a large range between the inner and outer disc radii, of at 
least -Rout/fiin ~ 10 4 . For these reasons, we are only able to 
obtain solutions using the time dependent code for high val- 
ues of f3 (/3 = 1.5), when the viscous time is a weak function 
of increasing radius. 



3 STEADY-STATE DISC SHAPE 

3.1 Comparison with the analytic solution 

Figure ^ shows the numerical steady-state solution for the 
f3 = case, and compares it to the approximate ana- 
lytic solution. We adopt units in which Ri n = 1, and take 
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Figure 1. Comparison of the numerical solution to the full 
warped disc equation with the approximate solution derived by 
Scheuer & Feiler (1996). The parameters are j3 = 0, i^io = ^20 = 
1, and uip = 200. The disc at i?, ou t = 10 s Ri n is fixed at an in- 
clination angle of 0.1 with respect to the spin axis of the black 
hole. The upper panel compares the inclination as a function of 
radius in the exact solution (solid line) and the analytic solution 
(dashed line). The lower panel shows the same comparison for 
two components of the tilt vector, l x and l y . The agreement is at 
the level of a few percent. 



fro — ^20 = 1. We choose ui p — 200, which gives an in- 
ner aligned region extending out to around 10 2 R ln . This is 
roughly appropriate for plausible AGN parameters, though 
we defer until later consideration of detailed disc models. 
Convergence to the limiting inclination at infinity is slow, 
so it is necessary to impose the outer boundary condition 
of i — 0.1 at a large radius, R = 10 s Ri n - At this radius l x 
and l y are chosen to facilitate comparison with the analytic 
solution of equation (^). 

Evidently, the analytic solution provides an excellent 
description of the disc shape. The inner disc is aligned to the 
rotation axis of the hole out to close to 10 2 Ri n , and warps 
to around half the limiting inclination by 10 3 R% n . Moving 
inward from infinity, the warp is twisted by an angle of ~ n 
before being flattened effectively into the aligned plane. 



3.2 Aligned radius 

For an annulus in the disc, the local timescales for precession 
and for transmission of warp are given by, 



V 20 RP ' 



We expect that the disc will be aligned with the spin axis 
of the hole at radii where the precession timescale is much 
shorter than the timescale over which the disc can commu- 
nicate warp inwards. The characteristic radius of the aligned 
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Figure 2. Disc shape for various values of /3. The upper panel 
shows the run of disc inclination i for models with (3 = 0, 0.1, 
0.2, 0.3 and 0.4 (larger /3 models are warped to smaller R). All 
models have u) p = 200 and uio = V20 = 1; i.e. the viscous time 
for both surface density evolution and warp is being kept fixed 
at the inner edge. The lower panel shows the extent of the inner 
aligned region of the disc as a function of [3, where -R a ii gn is defined 
implicitly via i(i? a ii gn ) = <5, with S = 0.01. The dashed curve 
shows the predicted relation if i? a ij gn scales with the radius where 
the viscous timescale for the warp equals the local precession 
timescale. 
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Figure 3. Radial dependence of the torque exerted on the disc. 
The upper two panels show the inclination i(R) and surface den- 
sity for three disc models, (i) i/i = V2 = 1 (solid line), (ii) 
v\ = 1, f2 = 8 (short dashed line) and (iii) v\ = 1/2 = R 1 / 4 , 
all normalised to the same accretion and precession rate. The 
units are chosen such that t vi and t„ 2 are also kept fixed at R ln . 
The lower panel shows where the resulting torques arise in the 
disc. We plot RAJ X and RAJ y for each model, where the total 
torque is given by the integral J AJdR. The RAJ X components 
are labelled in the lowest panel of the figure. For these parameters 



the greatest contribution arises at radii of ~ 10 
surface density and torque units are arbitrary. 



lO^ R 1U . The 



region -R a n gn will then scale with the radius where t p 
which is given by, 



-Ralien — C] 
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where Ci is a constant expected to be of the order of unity 
which we will use later to fit the numerical results. Note that 
this expression assumes both that the Lense-Thirring torque 
is strong enough to ensure that the inner disc is aligned, and 
that f3 > —1, so that t p increases with radius more rapidly 
than i„ 2 . 

Figure |^ shows how the disc shape varies with (3 for a set 
of models with fixed uj p and v\ = V2 = 1. From the numer- 
ical solutions, we compute i? a ii gn as the radius where i first 
exceeds some small inclination S. With Ci = 0.165 there is 
excellent agreement between the computed -R a ii gn and the 
scaling given by equation (^), i.e. the simple timescale ar- 
gument given above suffices to give an accurate idea of the 
scaling of the radius out to which the hole can enforce inner 
disc alignment. 



4 TORQUE 

The torque exerted on the misaligned disc as a result of the 
Lense-Thirring precession is given by the integral, 



dJ 



5. x L „ , „ 
p „, 2ttR<1R. 
H 6 



(9) 



Of course an oppositely directed torque of the same mag- 
nitude is exerted on the black hole. Computing J for the 
disc model shown in Fig. |l| we find that the torque agrees 
with that obtained from the Scheuer & Feiler solution at 
the level of a few percent (the small discrepancy being pri- 
marily due to the ignored inner boundary condition in the 
analytic solution, which leads to an error in the surface den- 
sity in the warp region). The torque scales as (wpi^) 1 ^ 2 M , 
again as found analytically. We have computed additional 
models with large inclination angles at infinity (up to 80° 
angle of misalignment) and for those models there is a larger 
discrepancy between the numerical and analytic results, but 
still only at the tens of percent level. 



4.1 Radial dependence 

Figure ^ shows the tilt angle and surface density profile for 
three fiducial models, all normalised to the same accretion 
rate. The solid line is computed for the parameters used for 
Figure [j], namely /3 = and v\o = f2o- The short dashed 
line maintains /3 — 0, but has 1/10 = 1, V20 = 8. Natarajan & 
Pringle (1998) find that consideration of the hydrodynamics 
of a viscous disc suggests that in AGN, 

1 (10) 
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in which case this ratio of viscosities would be appropriate 
for a plausible Shakura-Sunyaev a of around 0.25. Finally 
the long dashed line illustrates the effect of a declining sur- 
face density profile, plotted is a model with (9 = 1/4 and 
v\ — U2 . All three of these models are chosen to have the 
same viscous timescale (both for the surface density and for 
warp) at the (arbitrary) choice of inner edge radius. In con- 
ventional thin disc models (Shakura & Sunyaev 1973), fixing 
the viscous time at a given radius corresponds to a fixed disc 
thickness H/R at that radius. We note that in such models 
the value of /3 is fixed by solving for the vertical disc struc- 
ture, i.e. it is not a 'parameter' that can be varied but rather 
a property of the disc. In Section 5 we discuss the appropri- 
ate value for (5 in the region of the disc where the dominant 
torque arises. Also shown for each of the models is where 
radially the greatest contribution to the torque arises, we 
plot the components of the integrand in equation ([]) (the 
extra factor of R in the plotted quantity takes account of the 
logarithmic radial scale to give a true impression of where 
the strongest torque on the disc is found). 

For all the models, strong torques arise from a broad 
region that extends from just outside the alignment radius 
out to one or two orders of magnitude larger radius. The 
large radial range used in these calculations is thus essen- 
tial to capture the torque accurately. Increasing the ratio 
of i>2 relative to v\ reduces the timescale for diffusion of 
warp, t V2 = _R 2 /^2, and allows the warp to push in closer 
to the black hole. At fixed accretion rate increasing V2 even 
by this large factor leaves the surface density profile almost 
unaltered, and thus the net effect is simply to increase the 
magnitude of the torque and move it to smaller disc radius. 
Larger values of j3 likewise decrease the warp timescale at 
large radius and lead to a shrinkage of the aligned region. 
For the parameters used here, the disc shape is in fact rather 
similar for the two models in which U2 > 1 at large radius. 
However, for a given accretion rate, models with higher f3 
have lower surface density at large radii, which tends to de- 
crease L in equation (^ and reduce the integrated torque. 
Both these effects have a significant impact on the calcula- 
tion of the total torque. 



4.2 Scaling with f3 

The numerical results show that the Scheuer & Feiler (1996) 
solution, equation (^), provides an accurate description of 
the disc shape and the resultant torque on the black hole. 
For this solution, the torque integral in equation Q has a 
magnitude given by, 



V2KM , 
3^io 



(GM) 1/2 (u p v 20 
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(11) 



Guided by the numerical results, a straightforward extension 
of this result to non-zero /3 is possible. We first assume that 
the radius R w out to which we need to integrate equation 
(H) to find the torque scales as does the aligned radius, 
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If we additionally assume that the disc shape is primarily a 
function of R w , then the form of the integral in equation (^|) 
suggests the following scaling with f3, 
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Figure 4. Scaling of the torque with /3 at fixed mass accretion 
rate, where the viscosity v oc i?' 3 . The points show numerical 
results, the curves the fitting formula described in the text. The 
main parameters are: (i) solid line uio = V20 = 1, = 200, (ii) 
dotted line viq = U20 = 0.1, ui p = 200, (iii) short dashed line 
^10 = 1, V20 = 2, Up = 200, (iv) long dashed line viq = V20 = li 
u! p = 125. The torque on the hole drops by around an order of 
magnitude over the range of (3 considered here. 
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which reduces to the formula in equation ( Jll| ) if y3 = 0. We 
have allowed ourselves a single adjustable parameter in this 
expression, to best match the numerical results we adopt 
C 2 = 0.55. 

Figure ^ shows how equation (^) compares to the nu- 
merical results for choices of parameters testing the scalings 
implied by the formula. Excellent agreement (at better than 
the percent level) is obtained for two models with v\ = 7/2 
and varying oj p , while a single model computed with the 
time-dependent code at /3 = 1.5 also agrees well with the 
fitting formula. A model with V2 > V\ displays some sys- 
tematic error, caused by the influence of the increased ratio 
of V2/V1 on the surface density profile in the warp region, 
but equation ( p^ ) still provides a reasonable approximation 
to the numerical values for the torque. The main result is 
that the torque for a fixed accretion rate (and fixed viscous 
time at R ln ) drops by roughly an order of magnitude for disc 
models with f3 varying between (3 = and (5 = 1.5, while 
the dependence on other parameters (ui p , V20) remains close 
to that found in the (3 = case. 



5 ALIGNMENT TIMESCALE 

The torque on the black hole due to the interaction with the 
warped disc causes it to precess and become aligned with 
the outer disc's angular momentum vector on a timescale 
that, to factors of order unity, is just t a ii gn = |J|A'/| J| (e.g. 
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Figure 5. The timescale for the spin of a 10 s Mq black hole 
to become aligned with the angular momentum of the outer disc, 
plotted as a function of the mass accretion rate in units of the Ed- 
dington rate. The dashed curve is for V2 = v\, the solid curves are 
for v^/vi = l/(2o 2 ) with (from top downwards) a = 0.3,0.2,0.1. 

Scheuer & Feiler 1996). |J| is specified for a given disc model 
by equation ([l3|), while the angular momentum of a hole 
of mass M and spin parameter a is |J| = acM(GM/c 2 ). 
The resulting expression for t a ii gn is similar, but not equiv- 
alent, to the estimate given by Natarajan & Pringle (1998), 
taiign = (jJ|/|Jdisk|)(-Rw a r P / w p)> where Jdisk is the angular 
momentum of the disc within a warp radius 7? W ar P defined 
similarly to equation ([[2]). Our expression takes more com- 
plete account of the need to integrate the torque for all radii 
less than R w . For models of the disc where most of the an- 
gular momentum within the warp radius lies close to that 
radius, the final inferred timescale for alignment is similar. 

Obtaining an estimate for /3 and z/io requires a model 
for the vertical structure of the disc, which serves the pur- 
pose of relating the unknown central disc temperature, that 
controls the viscosity, to the effective temperature fixed by 
the requirement that viscous dissipation balance radiative 
losses. At the radii of the warp, gas pressure and an elec- 
tron scattering opacity k = 0.4 cm 2 g _1 are dominant. In 
this regime, Collin- Souffrin & Dumont (1990) obtain a ra- 
dial profile of column density N25 (in units of 10 25 cm~ 2 ) 
described by, 




where e is the radiative efficiency, L — eMc 2 is the bolomet- 
ric luminosity, I/Bdd = 47rc 3 i? 9 /ft is the Eddington luminos- 
ity, and R g = 2GM/c 2 . For this model, f3 = 0.6. Although 
this disc model includes the effects of irradiation of the up- 
per layers of the disc by the central source, we note that 



the scaling E oc a~ 4/5 M 3/5 M 1/5 i?~ 3/5 is identical to that 
of non-irradiated models (e.g. Sincell & Krolik 1998). The 
normalisation is within a factor of two. Other uncertainties 
thus dominate over that arising from different treatments of 
the vertical disc structure. 

Figure ^ shows the calculated alignment timescale for 
a 10 s Mq black hole accreting at rates between 10 -2 MEdd 
and MEdd, where MEdd = 4ircR g /(eK). We take a radiative 
efficiency e = 0.1, and assume a maximal Kerr hole, a = 
1. For these parameters, R w as defined in equation (jl2|) is 
Rw ~ 10 3 R g for reasonable ratios of U2/V1, so the main 
warped region of the disc falls self-consistently within the 
regime where gas pressure and electron scattering opacity 
are dominant. The Figure shows results calculated assuming 
that V2 = Ui, along with the expectation if V2/V1 = l/(2a 2 ) 
for a range of plausible a. The latter models of course predict 
more rapid realignment, as noted previously (Natarajan & 
Pringle 1998). 

For holes accreting at rates of the order of the Edding- 
ton limit, the alignment timescale is found to be 10 5 - 10 6 
yr. This is short compared to lifetimes of AGN, irrespective 
of what one assumes for the ratio of V2/V1. Producing the 
giant structures observed in radio galaxies may require an 
active epoch of the order of 10 8 yr, while some estimates for 
the quasar lifetime are similar (e.g. Haehnelt & Rees 1993, 
though these estimates are subject to considerable uncer- 
tainty). In these systems we expect that there was ample 
time for a rotating black hole, whatever its initial spin axis, 
to become aligned with the angular momentum of the disc. 
Thus jets accelerated from the inner disc should be perpen- 
dicular to the plane of the outer disc, and constancy in the 
direction of such jets implies a corresponding stability in the 
disc plane over time. 

Figure 6 plots the alignment timescale as a function of 
black hole mass, for a range of accretion rates. The align- 
ment timescale is found to be only weakly dependent on the 
mass of the hole, as expected given the M _1//16 dependence 
derived by Natarajan & Pringle (1998). For this figure we 
have assumed that v\ = v%, and that a = 0.1. Shorter align- 
ment timescales are of course predicted if V2 > as shown 
in Figure 5. 

For holes accreting at lower rates (relative to the Ed- 
dington limit), the timescale for alignment grows. It is not 
unambiguously shorter than estimates for the active phase 
of a hole in a galactic nucleus, and the answer depends cru- 
cially on differences in the disc's response to shear and warp. 
It remains true that i a ii g n <§C igrow, where i gr ow = M/M, for 
any M, so that alignment will always occur if an accretion 
event lasts long enough to contribute significantly to the 
hole mass. Perhaps more probably, however, low luminos- 
ity AGN may accrete negligible fractions of the hole mass 
in many brief episodes. In this scenario, it remains possi- 
ble that the hole spin accumulated over the whole accretion 
history will be sufficient to control the jet direction during 
subsequent low luminosity outbursts. The outer disc would 
then not be expected to be generally perpendicular to the 
jet direction. 
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Figure 6. The timescale for the spin of the black hole to become 
aligned with that of the outer disc is plotted as a function of 
the mass M of the black hole. From top downwards, the curves 
correspond to M = 10~ 2 M Edd , M = 10 _1 M B dd, and M = M Edd 
respectively. We have assumed a maximally rotating hole, an a = 
0.1, and V2 = v\- 



6 DISCUSSION 

In this paper we have considered the interaction between a 
misaligned accretion disc and a rotating black hole for the 
parameters appropriate to Active Galactic Nuclei. A com- 
bination of forced precession and disc viscosity allows the 
black hole to force the inner disc into a plane perpendicular 
to the hole's spin axis, but the resulting warped disc ex- 
erts a torque that eventually aligns the rotation of the hole 
with the angular momentum of the outer disc. We have pre- 
sented steady-state numerical solutions to the full warped 
disc equations, verified that the analytic solution of Scheuer 
& Feiler (1996) provides an accurate representation of the 
disc shape and torque on the hole, and generalised the re- 
sults for a range of disc models with varying viscosity. The 
range in models we consider lead to almost an order of mag- 
nitude variation in the integrated torque on the hole. 

We have computed the timescale for alignment to oc- 
cur, and find that for holes accreting at rates of the or- 
der of the Eddington limit (for a 10 8 Mq hole, this implies 
-^boiomotric = 2.5 x 10 46 erg/s, so these are luminous AGN) 
the timescale is short, 10 5 - 10 6 yr. This is the case for a 
range of assumptions as to the ratio 1*2/1/1, which character- 
izes the disc's relative response to azimuthal shear and to 
warp. The short derived timescale implies that the spin of 
black holes in such AGN rapidly adjusts to match the angu- 
lar momentum of accreting gas, even if the hole gains only a 
small (<C 10%) of its mass from accretion. For low luminos- 
ity AGN the prediction is more model-dependent, but the 
timescale for alignment in these systems is sufficiently long 
that the spin axis of the hole may be able to remain stable 
if individual accretion episodes are brief. Here the assumed 
value for vijv\ is crucial, and efforts to determine this ratio 



for specific angular momentum transport mechanisms and 
disc models would be valuable. 

Some powerful radio galaxies display jets that appear 
to have maintained their direction for long periods of time. 
Our results imply that this cannot be due to any intrinsic 
stability imparted by the spin of the black hole, but instead 
must reflect a long term constancy in the angular momentum 
of the outer regions of the accretion disc. If the gas derives 
from a reservoir set up by a single accretion event then this 
constancy would not be surprising. Alternatively, a preferred 
axis for gas arriving in the nuclear regions might be the 
consequence of interactions between inflowing gas and the 
galactic potential. 

We have calculated the torque assuming that the disc 
is able to relax into and maintain a steady-state configura- 
tion as realignment occurs. This assumption is justified in 
the strongly warped inner region, where t vi and t V2 are in- 
deed smaller than the derived i a iign, but fails at large radii 
of 10 4 - 10 5 R g where our solutions show that the disc main- 
tains a small but significant warp. This will not increase the 
alignment timescale (before the disc reaches a steady-state 
the torque in our time dependent calculations is larger than 
the limiting values we solve for), but does suggest that fol- 
lowing alignment of the hole the outer disc might retain a 
modest warp for a long period. For a 10 8 Mq hole a radius of 
10 4 R g corresponds to ~ 0.1 pc, a scale which is observable 
with VLBI in nearby AGN such as NGC4258 (eg. Miyoshi et 
al. 1995). Some of the warps in maser discs could then sim- 
ply reflect the slow decay of the initial conditions of a past 
accretion event, and be devoid of a persistent forcing mech- 
anism. Of course NGC4258 itself has an accretion rate much 
lower than the models considered here (Gammie, Narayan 
& Blandford 1999), and self-gravity is likely to be important 
in its disc (Papaloizou, Terquem & Lin 1998), but the main 
point — that t„ 2 at a few tenths of a parsec is likely to be 
long — remains true. 

A related question is whether the disc at R ~ 10 2 R Q 
is able to attain a steady state, even in the absence of a 
warp. This region of the disc may be thermally unstable 
(Lin & Shields 1986; Clarke 1989; Siemiginowska, Czerny & 
Kostyunin 1996) and prone to a limit cycle involving large 
changes in the accretion rate. The interplay between such a 
cycle and the dynamics of a warped disc would be complex 
and intriguing, though a more promising setting for explor- 
ing such a scenario would be in galactic black hole systems, 
if those possess warped discs. 

For AGN with rapidly spinning holes, the main warped 
region occurs at radii of ~ 10 2 R g . During realignment 
episodes there is a substantial deposition of energy, origi- 
nating in the spin energy of the hole, into this region of the 
disc. This is particularly the case if V2 vi, as we have 
argued is likely in this paper and previously (Natarajan & 
Pringle 1998). The disc in the vertical structure models we 
have considered is quite thin at these radii and can easily 
radiate locally a significantly enhanced luminosity, but sub- 
stantial changes to the vertical structure might be expected 
(note that we have assumed that the vertical structure of 
a warped disc is just like that of a planar disc, which is 
wrong on several counts). Moreover, the enhanced luminos- 
ity would increase the disc flux at optical and infrared wave- 
lengths, which would in any event be raised due solely to the 
greater covering fraction of the warped disc as seen from the 
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central source. Since the alignment timescale is short, sys- 
tems of this kind would be rare — most AGN would harbour 
quite flat central discs — but potentially bright. 
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APPENDIX 

The ODEs for the steady state shape of the disc are most 
easily solved by expressing equation (Q) in terms of separate 
equations for the tilt vector I and surface density E. For a 
Keplerian potential these are (Pringle 1992), 
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where the primes denote derivatives with respect to R, and 
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Setting the time derivatives to zero these constitute 3 second 
order equations for E and any two components of /. We solve 
these with the boundary conditions described in Section 2.3 
using the NAG routines D02GAF and D02RAF, which im- 
plement the finite difference technique with deferred correc- 
tion described by Pereyra (1979) (for a general introduction 
to such methods, see e.g. Press et al. 1992). An iterative 
approach is employed, in which we first solve equation ( |l6| ) 
assuming zero tilt, and then use the resulting surface den- 
sity profile in the solution of equation ( |j"i|) . The solution for 
the tilt vector / is then recycled for use in equation (|l6|), 
and we loop until convergence is achieved. Typically only a 
small number of iterations are required, since the changes 
to the surface density profile for even quite strong warps 
vary smoothly with changing disc shape. A finely space fi- 
nite difference mesh is required to obtain a solution using 
this scheme, we use between 8000 and 16,000 mesh points 
evenly spaced in log R. Even so, we have found that ob- 
taining solutions with this scheme is still rather difficult, 
especially when /3 > or v\ / U2- For these cases, we start 
with the easy /3 = 0, v\ = i>2 solution, and step towards the 
desired parameters using the previous solution at each step 
as the initial guess. 

We have been unable to obtain solutions using the above 
method for large values of /3. In this regime, we instead 
evolve directly equation (|l|) using the numerical code de- 
scribed by Pringle (1992) until a steady-state is obtained. 
This code uses an explicit first order finite difference tech- 
nique, which conserves angular momentum to machine pre- 
cision. We have modified the boundary conditions to be as 
close as possible to those described in Section 2.3, so that 
the results are directly comparable between the two solution 
methods. The time dependent code is vastly more expensive 
to run, for the runs described here we are restricted to 100 
radial grid points, again spaced evenly in log R. Neverthe- 
less this still provides reasonable accuracy for computing the 
torque. 



© 0000 RAS, MNRAS 000, 000-000 



